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Abstract 

We propose a new fractional statistics for arbitrary dimensions, based on 
an extension of Pauli's exclusion principle, to allow for finite multi-occupancies 
of a single quantum state. By explicitly constructing the many-body Hilbert 
space, we obtain a new algebra of operators and a new thermodynamics. The 
new statistics is different from fractional exclusion statistics; and in a certain 
limit, it reduces to the case of parafermi statistics. 
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The thermodynamics of a macroscopic system is determined microscopically by 
the statistics of its constituent particles and elementary excitations. Herein lies a 
fundamental significance of statistics. Ever since Heisenberg's second paper on matrix 
mechanics, it has been known that a many-body wavefunction is symmetric under 
permutations of identical bosons, but it is antisymmetric for identical fermions. The 
corresponding commutation and anticommutation relations bilinear in field operators 
result in bose and fermi statistics respectively. Particles are accordingly classified 
into bosons and fermions. The overriding difference between the two groups is that 
bosons condense while fermions exclude. But it is natural to inquire whether there 
are any meaningful generalizations of statistics intermediate between these two. 

Attempts to generalize statistics dates back at least to Green's work in 1953 
I] || . Green found that the principles of quantum mechanics also allow two kinds of 
statistics called parabose statistics and parafermi statistics of positive integral order 
M (the M=l cases reduce to the familiar Bose-Einstein statistics and Fermi-Dirac 
statistics respectively). They are described by trilinear commutation relations among 
the creation and annihilation operators. Subsequently, the case of non-integral M was 
investigated for possible deviations from Bose and Fermi statistics, and in particular, 
for possible violations of Pauli's exclusion principle 0][|J. This saga culminated with 
a recent study of infinite statistics |JJ in which all representations of the symmetric 
group can occur; this case is realized by the q-mutator algebras. 

Another type of interpolating statistics, spearheaded by Wilczek, is provided by 
the concept of anyons ||. Anyons are particles whose wavefunctions acquire an ar- 
bitrary phase when two of them are braided; they obey fractional statistics. More 
recently, Haldane introduced another definition of statistics based on a generalization 
of the Pauli principle J7j||. Unlike the anyon fractional exchange statistics which is 
meaningful only in two spatial dimensions, Haldane's fractional exclusion statistics is 
formulated in arbitrary dimensions. The thermodynamics based on exclusion statis- 
tics is studied in Ref. ||. The issue whether anyons obey fractional exclusion statistics 
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in the framework of quantum field theory is addressed in Ref. 0. In this Letter we 
introduce another type of statistics which is based on a (different) logical extension 
of the Pauli principle. In the new statistics, we allow for multi-occupancy of a single 
quantum state by up to a maximum number M of identical particles. The M = 1 case 
yields the conventional Fermi statistics while the M — > oo case corresponds to the 
conventional Bose statistics. The new statistics bears some resemblance to fractional 
exclusion statistics but is distinct from it. 

In the literature a frequently used approach to study the quantum features of 
(especially fermi and fermi-like) many-body systems is to start with assumptions of 
commutation relations among certain operators. In this letter, we would like to start 
with the construction of Hilbert space of quantum states and then "derive" some 
relations of the operators. Consider the Hilbert space spanned by the eigenstates of 
the particle number operator N, {\j >, j = 0,1,2,- • ■}. The j-particle state \j > 
satisfies 

N\j>=j\j>, <j\k>=5 jk . (1) 

Obviously, the ground state is the zero-particle state |0 >. We define the one-particle 
state as a superposition of M single states |e* >: 

|1 >= cic 2 • • • c M -i|ei > +sic 2 • • • c M _i|e 2 > H hs M -i|e M > , (2) 

where q = cos6*i, Si = sin^j, 9i ^ 0, or multiples of ir/2. The use of the angles is just 
a matter of convenience. In this case, the one-particle state can be thought of as a 
unit vector on a (M— l)-sphere. The M single states |ej >, i = 1, 2, • • • M, satisfy 

< 6i\ej >= 5ij , and (|e» >) 2 = . (3) 

The latter condition can be understood as a reflection of the Pauli exclusion principle, 
namely, a two-|e^ > state is forbidden. 

The two-particle state is defined as a superposition of the tensor products 
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M 

|2 >=J2c ij \e i > \ej > . (4) 

i<j 

Similarly, three and more particle states can be defined as a superposition of tensor 
products |e, > \tj > > with c- number coefficients and so on. How to determine 
Cij, c^k and so on, in term of the angles 0j, will be discussed later. We should 
emphasize that in a tensor product of single states, we adopt the rule that the order 
of single states does not make a difference, i.e. |e$ > \ej >= |e 3 - > |e, >. Due to the 
second condition in Eq.(|3|), the M-particle state 

M 

im >= n h > > ( 5 ) 

is the maximal-particle state, the state with the maximum number of M particles. 
There are no states beyond it, i.e. states \j > for j > M do not exist. 
Now, we introduce the annihilation and creation operators, a and a\ 

a\j >= fj\j - 1 > , a f |j >= f* +1 \j + 1 > • (6) 

In particular, a\0 >= and a^|0 >= |1 >. Hereafter we choose fi and all fj to be real 
for convenience (this does not affect the physics we will discuss in the single particle 
species case). One can readily check the following commutation relation of N with a 
and at 

N a - aN = -a , and Na ] - a ] N = a ] . (7) 
Starting from the one-particle state Eq.(^), using the relation 

( a t)J'|0 >= (|1 >) j ~ \j > , or a j \j >~ |0 > , (8) 

the normalization condition in Eq.(^j), and / = and j\ = 1, one can systematically 
determine the amplitude fi and the coefficients c^, and c^k, and so on, set by set. 
To be concrete, let us take M = 2 and 3 as examples. For M = 2, the one-particle 
state is |1 >= c|ei > +s|e2 >. Acting the creation operator a) on |1 >, one has 
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at\l >= (|1 >) 2 = 2cs|ei > |e 2 >= / 2 |2 >. Therefore / 2 = 2cos#sin6'. Obviously |2 > 
is the maximal-particle state since |3 >= 0. This is tantamount to fj = for j > 2, 
or equivalently (at) 3 = a 3 = 0. 

For M = 3, the one-particle state is |1 >= cic 2 |ei > +sic 2 |e 2 > +s 2 |e 3 >• 
From at\l >= (|1 >) 2 = / 2 |2 >, we obtain f 2 = 2s 1 CiC 2 ,^l + t|/ s i c i with t 2 = 
tan6> 2 and |2 >= 1/ ^1 + t\j s\c\ (|ei > |e 2 > +t 2 /si|ei > |e 3 > +t 2 /ci|e 2 > |e 3 >). 
Using at|2 >= / 3 |3 > with |3 >= \e\ > |e 2 > |e 3 >, we readily obtain / 3 = 
3s 2 /y / l + t|/ s i c i- The amplitudes /j = for j > 3, equivalently (a"!") 4 = a 4 = 0. 

Let us consider the operator algebra of a and a\ First, for a general M, there are 
no simple bi-linear operator relations like [a, a if ]± = 1 that the conventional boson 
and fermion operators satisfy. Instead, here we have 

(aat + ala)\j>=(ff + ff +1 )\j> . (9) 

For M = 1, we recover the anti-commutation relation for fermions. 

For any given M, a J = (a^)- 7 = only if j > M. Therefore, there exist (M+l)-linear 
relations between a and a*. For example, for M = 2, one has the cubic relations 

aV + / 2 W = / 2 2 a, (10) 
a 2 o) + a) a 2 + aa) a = Tr(aa^)a , (11) 

where Tr(aa^) = Y^f=\ fj, P ms the hermitian conjugate relations. 
For M = 3, there are the quartic relations like 

aV + / 3 2 ata 3 = / 3 V , (12) 
a 3 a^ + a j a 3 + a 2 a^a + aa^a 2 = Tr(aat)a 2 , (13) 

plus the hermitian conjugate relations. We speculate that the multi-linear relations 
among a and at in the limit M — > oo (and with a suitable choice of the theta angles) 
actually reduce to the bilinear commutation relation for Bose statistics. 



5 



Next we consider the particle number operator N in terms of the creation and 
annihilation operators a) and a. One way to do this is to ssume, for a given M, 

N = data + C 2 (a^) 2 a 2 + ■■■ + C M (a f ) M a M . (14) 

Then the M coefficients Cj, j = 1,2,- ■ -,M, can be determined by using the M 
independent equations N\k >= k\k >, k = 1,2, • • -,M. For example, for M = 2: 
C x = 1 and C 2 = (2 - fl)/fl For M = 3: C x = 1, C 2 = (2 - /f)//f, and 
C 3 = (3//| - 3 + /f)//f. 

The particle number operator iV can be expressed in other forms for certain values 
of the theta angles. For example, if the O^s in Eq.(Q) are chosen so that the one 
particle state is |1 >= l/\/M(|ei > + \e 2 > + ■ ■ ■ \eu >) and if the operators a and 
a) are replaced by operators 6/a/M and tf/y/M, the particle number operator takes 
the form N = \{b^b — btf) + yl. This latter form was used in the study of parafermi 
statistics 00. 

With the Hilbert space of quantum states for the new statistics and the diagonal 
particle number operator now available, it is natural and straightforward to con- 
sider the quantum statistical mechanics of a system that is compatible with such a 
construction (note that systems of this kind are not necessarily described by a free 
theory). 

Let us assume a single particle to have energy e = e(p), then the Hamiltonian 
operator takes the form 

H = eN . (15) 

The energy spectrum of the system for a given M is given by {e, 2e, • ••, Me}, similar 
to that for a spin system in a magnetic field. 
The grand partition function is 

M M 

Z = Tte-Ptf-vA) = J2< 3\e- p{R -» N) \] >= nE(^ e ) j > ( 16 ) 

j=0 P j=0 

where ft is the reciprocal of temperature T, \i is the chemical potential, the fugacity 
is z = e^, and the Boltzmann constant is — 1- For M = 1, Z — Y\ p (l + ze~@ e ) 
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is the partition function for free fermions; while for M = oo, Z = ]J p 1/(1 — ze~^~) 
recovers the partition function of free bosons [Kj. For M^l and oo, Z describes a 
system interpolating between free fermions and free bosons. Note that the resulting 
thermodynamics is insensitive to the particular values of the set of amplitudes fj. In 
particular, for the M = 2 case, even when the one-particle state |1 > is predominantly 
| ei > or |e 2 >, the statistics is very different from Fermi statistics, suggesting that 
the models discussed in Ref.[Q do not yield weak violations of the Pauli principle as 
correctly pointed out in Ref.[|J. 

With the grand partition function Z, one can calculate various thermodynamical 
quantities. The particle number 

Accordingly, the average occupation numbers are 

" (<) = tu^y ■ (18) 

At T = 0, n(e) = for e > //; while n(e) = M for e < //. The fermi energy ep is 
defined by the particle density n = N/V = (1/V) S P < PF n ( e ) a ^ absolute zero. As 
T -> oo, n(e) = M/2. 

One can also calculate the entropy S by applying S = —OF/ dT, where F = — TlnZ 
is the grand potential. Using Eq . ( 18 ) to invert ze~@ £ in term of the average occupation 
number n(e) one can then express S in term of n(e). For example, for M = 2, we find 

S = J2 (-lnn(e) + (1 - n{e))lnx + ln(l + 2x)) , (19) 



where x = (yl + 6n(e) - 3n 2 (e) - 1 + n(e))/2(2 - n(e)). 
The equation of state is given by 

M 

PPV = hxZ = J2 \nJ2(ze-?y , (20) 

P 3=0 

where P denotes the pressure and V the volume. In the large volume (V — ■> oo) limit, 
we replace the sum over momentum p by the integral over p: J2p ~ * V I d D p/ (27r) D . 



Such a replacement is clearly valid only if the summand is finite for all p. For the 
bose gas (the limit of M — > oo), the summand — ln(l — ze~@ e ) in Eq. (p0|) diverges as 
the fugacity z — *■ 1, because the single term corresponding to p = diverges. This 
is of course related to the Bose-Einstein condensation. On the other hand, for any 
finite M, the summand in Eq.(|20|) is finite for any value of e(p). 

In all our discussions so far we have made no reference to any specific spatial 
dimensions. We now consider a planar system. Furthermore, we assume the single 
particle energy e take the form e(p) = p 2 /(2m), with m being the (effective) mass 
of the particles or excitations so that the system is an ideal gas. Using Eq . fl20|) and 
Eq . (|T7|) , and performing the integrations over p, we readily obtain 

i oo k / 7 Mfc \ 

PP = C^ 1 ). ( 21 ) 



A 2 k 2 \ M + 1 
1 

Y 2 , „ . „ 

N 1,1- ^ (M+1) 
V = A^ ln ^7 



7T 2 M M„ N2 ™z~ k ( z~ Mk * 



(z>1), (22) 



n = 77 = ^— — , (23) 



where A = y27r/5/m is the thermal wavelength. Solving for z in Eq.(p3|) and substi- 
tuting it into Eq.(|21"D, in the high temperature and low density limit, i.e. \ 2 n 1, 
we can conduct a virial expansion in the form f3P = n (1 + B2\ 2 n + B 3 (X 2 n) 2 + •••). 
For M = 2, we find B 2 = -1/4, B 3 = 25/36, • • •; and for M = 3, B 2 = -1/4, 
B 3 = 1/36, ■ ■ ■. Actually from Eq. (pl|) and Eq . (p3|) , it is not difficult to check that 
for any M > 1, the second virial coefficient is —1/4, the same as for the ideal bose 
gas (M = oo). It implies that (for sufficiently small A 2 n) the quantum effect on the 
ideal gas, for all M except M = 1, is equivalent to an attractive "interaction" among 
excitations. For ideal fermi gas (M = 1) this effect is a repulsive one as B 2 = +1/4. 

In the low temperature and high density limit, i.e. X 2 n 3> 1, most particles are 
in the states with energy e < ep = 2im/M.m. Using U = F + TS + fiN, we obtain 
the internal energy 

"=i" e '( i+ sd£u£> ,+ -)- (24) 



The first term is the ground state energy, a result that can be verified by using 
J2p< PF Mp 2 /(2m) = Nepjl. From the above equation the specific heat at constant 
volume can be readily found: Cy/N ~ 3(m+i) These suggest that in the low 
temperature and high density limit, a system in which each quantum state has a 
maximum multi-occupancy of M < oo is like the fermion system (M= 1). 

A comparison of the quantum statistical mechanics obtained here with that for 
exclusion statistics can now be made. Both fractional statistics, based on general- 
izations of Pauli's exclusion principle, are well defined in arbitrary dimensions. At 
zero temperature, the n(e) distributions in the two thermodynamics are the same if 
the statistical parameter g in exclusion statistics is identified with 1/M in the new 
statistics. But there are profound differences between the two statistics. The second 
virial coefficient for a free planar exclusion statistical system, such as the free anyon 
system, is given by B2 = 1/4 — g/2 (where g = 0, 1 for fermion and boson statistics 
respectively) |§||. Accordingly the statistical interaction is attractive, neutral, or 
repulsive, depending on the value of g. In particular, for semions with g = 1/2, it 
is neutral (to this order); and for others with fractional g = 1/M (M > 2), it is 
repulsive. In contrast, we find an attractive interaction for all M > 1 in the new 
statistics discussed above. Moreover, the particle distribution n(e) are in general 
very different in the two statistics. For example, the distribution of semions from 
the exclusion statistical derivation takes the form ||: n(e) = 1/^1/4 + e 2 P e /z 2 . This 
is different (except at zero temperature) from the one given by Eq . fll8|) for M = 2, 
n(e) = (2 + e^ e /z)/(l + e^ € /z + e 2/3e /z 2 ). Integrating out the distributions over e, we 
find the resulting densities are different too. The statistical weights W for fractional 
exclusion statistics are also quite different from those for the new statistics. For in- 
stance, for N identical semions occupying G states, it reads W = 
in exclusion statistics || , whereas for M = 2 in the new statistics we obtain 



g + Un-i) 



N 
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where [f ] denotes N/2 and (JV-1) /2 for JV = even and odd, respectively. We conclude 
that the statistics associated with multi-occupancy of a single quantum state and the 
resulting operator algebra are different from the exclusion statistics defined in Ref.0 
(in contrast to a recent proposal |IT|| , in which the exclusion statistical parameter g is 



assumed to be connected to the maximum occupancy number M by g = 1 /M and the 
statistical distributions in exclusion statistics are connected to the amplitudes fj). 

We thank K. Dy E. Merzbacher, V.P. Nair, GW. Semenoff, and Y.-S. Wu for 
useful conversations. This work was supported in part by the U.S. DOE grant No. 
DE-FG05-85ER-40219. 

Note added: After this paper was completed, we noticed a recent interesting work 



12|, in which an issue relevant to this paper was also addressed. 
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